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Coding signals 



This invention relates to method of coding signals and to apparatus for storing, 
transmitting, receiving or reproducing signals. 



5 A common method of storing audio signals is to use parametric coding to 

represent audio signals, especially at very low bit rates, typically in the region from 6 kbps to 
90 kbps. Examples of the use of parametric coding used in this way are included in "Low bit 
J" rate high quality audio coding with combined harmonic and wavelet representation" in 
Q Proceedings of the IEEE International Conference on Acoustics, Speech and Signal 
fP Processing, Volume 2, pp 1045 to 1048, 1996; "Advances in Parametric Audio Coding" in 
I J Proceedings of the 1999 IEEE Workshop on Applications of Signal Processing to Audio and 
fy Acoustics, pp W99-1-W99-4, 1999; and "A 6 kbps to 85 kbps scalable audio coder" in 
1^ Proceedings of the IEEE International Conference on Acoustics, Speech and Signal 
!■« Processing, Volume II, pp 877-880, 2000. In these examples, a parametric audio coder is 
|| described, in which an audio signal is represented by a model, with parameters of the model 
- being estimated and encoded. These examples use a parametric representation of an audio 
signal based on decomposition of an original signal into three components: a transient 
component, a tonal (sinusoidal) component, and a noise component. Each component is 
represented by a corresponding set of parameters, as described in the three documents above. 
20 A transient component of an audio signal can be characterized as an isolated element of the 
audio signal which is relatively short lived, and is represented by a sharp increase in energy 
of the audio signal. 

It has been found that having a dedicated model for the transient component of 
an audio signal proves to be beneficial for parts of audio signals with sharp attacks, because 
25 sinusoidal and noise models cannot easily represent such perceptually important events and 
poor modeling can result in audible artifacts such as a pre-echo. A pre-echo occurs when the 
modeling error distributes the transient event to the samples before the transient beginning 
and when the resulted distortion is large enough to become audible. The distribution of the 
modeling error to the samples before the transient begiiming results from the segment-by- 
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segment analysis of an input signal in an audio coder. If a transient occurs in the middle of 
an analysis segment, then either a lot of coding resources are required in order to accurately 
model the transient, or the modeling error distributes to the whole analysis segment. 
Modeling error of the samples preceding a transient is typically perceptually more apparent 
5 than at samples after the transient, because of a weaker masking from the transient event 
itself 

In "Residual modeling in music analysis-synthesis" from Proceedings of the 
IEEE International Conference on Acoustics, Speech and Signal Processing, Volume 2, pp 
1005-1008, 1996 it is shown that transient components cannot satisfactorily be represented 

10 by sinusoidal and noise models alone. 

It has been shown previously in "Robust exponential modeling of audio 
signals" from Proceedings of the IEEE International Conference on Acoustics, Speech and 
Signal Processing, Volume 6, pp 3581-3584, 1998, that transients can be modeled efficiently 
C3 using sinusoids with exponentially modulated amplitudes (referred to below as damped 

1 1 sinusoids). In the text below damping coefficients can be any real number, and positive 
j^^ values correspond to increasing amplitudes rather than to truly decreasing amplitudes. In 

fU "Robust exponential modeling of audio signals" (see above) an audio signal was analyzed on 

a segment-by-segment basis and each segment was represented as a sum of damped 
^* sinusoids. A problem arises with this type of coding when a transient starts in the middle of a 

a 

2Q given segment. Compared to the case where transient starts in the beginning of a segment, 
^ the number of damped sinusoids needed to model the transient well increases considerably. 

If a transient is not modeled properly, the modeling error is distributed over the whole of a 

given segment resulting in audible pre-echoes. 

In the MPEG-1 Layer III audio coding algorithm, as described in "ISO- 
25 MPEG-1 Audio: a generic standard for coding of high-quality digital audio" in the Journal of 

the Audio Engineering Society, Volume 42, pp 780-792, October 1994. The segmentation is 

defined simply by the lengths of the long and short windows. 



30 It is an object of the present invention to address the above mentioned 

disadvantages. To this end the invention provides a method of coding and an apparatus for 
coding as defined in the independent claims. Advantageous embodiments are defined in the 
dependent claims. 
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According to a first aspect of the present invention the coding of an input 
signal comprises: 

- estimating a location of at least one transients in a time segment of the input signal; 

- modifying the location of the transient so that the or each transient occurs at a specified 
location on a predetermined time scale to obtain a modified signal; and 

- modeling the modified signal. 

The use of restricted time segmentation in the form of a specified location on a 
predetermined time scale to provide the only locations for the transients advantageously 
reduces the number of bits needed to describe the segmentation. Also the modification 
procedure has lower computational cost compared to a full precision segmentation procedure. 

Each transient is preferably re-located to a nearest specified location of a 
plurality of possible locations on the predetermined time scale. 

The specified locations on the predetermined time scale may be defined by 
integer multiples of a predetermined minimum time segment size. The predetermined 
minimum time segment size may have a length in the range of approximately 1 millisecond 
(ms) to approximately 9 ms, most preferably in the range of approximately 4 ms to 
approximately 6 ms. 

The use of a restricted time segmentation as described advantageously 
simplifies the modeling procedure significantly, if rate-distortion control is used to distribute 
coding resources between teansient, sinusoidal and noise components of the input signal 
being modeled. 

The modeling preferably uses damped sinusoids. 

The audio signal is preferably sampled at a rate of approximately 5 to 50 kHz, 
most preferably 8, 16, 32, 44.1 or 48 kHz. The video signal is preferably sampled at a rate 
of approximately 5 to 20 MHz. 

The restricted time segmentation may also be applied to tonal and/or noise 
components of an input signal. 

The estimation of the location of transients may be carried out using an 
energy-based approach, preferably with a moving window method, most preferably using 
two sliding windows. 

The use of an energy-based approach allows the advantageous estimation of 
both very short transients and longer transients. 
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The location of transients may involve the location of a beginning said an end 
of each transient. 

Preferably each located transient is moved by a cut and paste method from its 
original location to begin at a location on the predetermined time scale. 

The cut and paste method simply removes that part of the input signal 
identified as a transient and moves it to the new location. Thus the step is very simple to 
implement. 

A remaining section of the input signal between two located and modified 
transients is preferably time- warped to fill the gap remaining following the relocation. The 
time-warp may be a lengthening or a shortening of said remaining section. 

By using knowledge of soimd perception, including pitch perception and 
temporal masking effects, the time-warping is a simple method with which to restore the 
remaining signal after modification of the transients. 

The time-warping preferably preserves the amplitudes of edge-points of the 
modified signal, preferably by a band limited interpolation method. 

The time-warp is preferably carried out by interpolation where the change in 
the fundamental fi-equency,/o, of the remaining section is less than approximately 0.3%, most 
preferably less than approximately 0.2%. 

Otherwise, the remaining section is preferably split in to a first length 
immediately after the modified transient and a second length. Preferably, the first length is 
approximately 8 ms to 12 ms, most preferably approximately 10 ms. The first length is 
preferably interpolated if the change of fundamental frequency caused is no more than 
approximately 1 .6% to 2.4%, most preferably no more than approximately 2%. For the 
second length, the change of fundamental frequency is preferably not more than about 0.16% 
to 0.24%, most preferably approximately 0.2%. 

Where the interpolation is insufficient to fill a gap in the remaining section an 
overlap-add procedure is preferably used. 

The modification of the location of the or each transient may be performed 
using a transformation into a frequency domain, preferably with a discrete cosine transform. 
The resulting sinusoidal representation may then be analyzed for transient locations using a 
Harming window. Preferably, the Hanning window has a length of approximately 512 
samples (where a sample has a length of one divided by a sampling frequency of the input 
signal), preferably with an overlap between Hanning windows of 256 samples. 
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The input signal is preferably processed by dividing the input signal into a 

plurality of time segments. The time segments may have a length in the range of 

approximately 0.5 s to 2 s, preferably a length of approximately 1 s. 

Adjacent time segments are preferably arranged to overlap, preferably by 
5 approximately 5% to approximately 1 5% of their length, more preferably the overlap is 

approximately 10% of the time segment length, which overlap may be approximately 0.1 s. 

Where a transient is located in an overlap of the adjacent time segments, the transient 

location is modified in the time segment in which the transient is most centrally located. 

The provision of an overlap in adjacent time segments advantageously allows 

10 the selection of the time segment in which the transient is most centrally located, or more 
importantly furthest from the beginning or end of the time segment. 

The invention extends to decoding audio or video signals coded according to 
f J the coding of the first aspect. 

An apparatus according to an embodiment of the invention may be an audio 

11 device, e.g. a solid state audio device. 

[fi AH of the features disclosed herein can be combined with any of the above 

^ ^ aspects, in any combination. 

N= Preferred embodiments of the invention of the invention provide coding 

y... 

signals which coding has a more simplified analysis procedure than has previously been 
M described, coding signals which coding has a lower computational cost than equivalent 
U methods, and coding signals which coding results in a reduction of the number of bits needed 

to describe a segmented signal. 

Additional side information may be included in the bitstream to dewarp the 

signal at the decoder side. With the appropriate dewarping, temporal misalignment of stereo 
25 signals can be avoided. 



Specific embodiments of the present invention will now be described, by way 
of example, and with reference to the accompanying drawings, in which: 

Figure 1 shows the performance of a damped sinusoidal model in the case of a 
restricted segmentation of an audio signal for an original and a time shifted transient for a 
first embodiment; 

Figure 2 shows an original transient and its reconstruction with 25 damped 

sinusoids; 
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Figure 3 shows a time shifted transient and its reconstruction with 25 damped 
sinusoids for the first embodiment; 

Figure 4 is a flow diagram of the steps involved in the method of coding audio 
signals in the first embodiment; 
5 Figure 5 is a diagrammatic illustration of the modification of transient location 

in a second embodiment; 

Figure 6 is a diagrammatic illustration similar to that of Figure 5; 
Figure 7 shows an original transient and its reconstruction; 
Figure 8 shows a shifted transient and its reconstruction according to the 
10 second embodiment; 

Figure 9 is a flow diagram of the steps involved in the second embodiment; 

and 

fn Figwe 10 is a schematic diagram of an audio encoder and an audio decoder 

j3 utilizing the methods described herein. 

I 

m 

ft^ The first method disclosed herein, and as shown in Figure 4, uses a restricted 

|,.^. time segmentation, in which segments of an audio signal are defined by integer multiples of a 
p predefined minimvim segment size, which in the example used is 5 ms, but of course this 
^ could vary. In view of the restricted time segmentation the transient component of the audio 
r~ signal is modified such that transients can start only at the beginning of a segment. The 

modified signal is then modeled, in this example by using damped sinusoids. This results in 
an efficient representation of transients with damped sinusoids. 

The coding of audio involves a first step of modifying the location of transient 
25 elements of the signal so that the transients can occur only at locations defined by a relatively 
co^se time grid, as described below in the discussion of experimental results. In order to 
modify the locations of transients in the audio signal the following steps are taken: 
1 . The transient component of an original audio signal is estimated and is subtracted from 
the original audio signal to form a residual signal. 
30 2. The locations of the estimated transients are then modified in such a way that the 
transients can only occur at locations specified on a grid. 

During the transient estimation and modification, it has been verified that 
when the modified transient signal is added to the residual signal obtained in step 1 above, 
there is no perceptual difference between the obtained signal and the original audio signal. 
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In order to modify the transient locations it is necessary to estimate the 
transient component of the original audio signal to be coded. It is possible to use different 
transient models in parametric coding of audio. One example which has been used is the 
transient model based on duality between the time and frequency domain presented in 
"Transient modeling synthesis: a flexible analysis/synthesis tool for transient signals", in 
Proceedings of the Intemational Computer Music Conference, pp 25-30, 1997. 

In more detail, the transient estimation model presented in the above reference 
is based on the duality between the time and the frequency domain. A delta impulse in the 
time domain corresponds to a sinusoid in the frequency domain. Furthermore, a sharp 
transient in the time domain corresponds to a frequency domain signal which can be 
represented efficiently by a sum of sinusoids. More specifically, the transients are estimated 
using the following steps. 

1 . A discrete cosine transform (DCT) is used to transform a time domain segment to the 
frequency domain. The segment size (equivalently, the DCT size) should be sufficiently 
large to ensure that a transient is a short event in time (thus, transformed to the frequency 
domain, it can be modeled efficiently by sinusoids). A block size of about 1 s has been 
found to be sufficient. 

2. The frequency domain (DCT domain) signal is analysed with a sinusoidal model. One 
example which has been used is a consistent iterative sinusoidal analysis/sjnithesis with 
Hanning-windowed sinusoids, as described in "High quality consistent analysis-synthesis 
in sinusoidal coding", from Proceedings of the Audio Engineering Society 17* 
Conference "High quality audio coding", pp 244-250, 1999. 

The sinusoidal analysis of a DCT domain segment is done on a segment by 
segment basis. As a result, the DCT-domain segment is represented as 



/ = 0,...,Z,-l,i=l,...,/ 

where L is the length of the sinusoidal segments (the shift between sinusoidal segments is 
L/2). The length of the sinusoidal segments, L, is a small fraction of the DCT size, N. h(l) 
are samples of the Hannmg window, and {Aij, co,j, 0ij} are amplitudes, frequencies and 
phases of the estimated sinusoids respectively. The index i denotes a particular sinusoidal 
segment within the DCT-domain segment, while the index j denotes a particular sinxisoid 
within the sinusoidal segment. The information about the location of a transient in a time 
domain segment is contained in the frequency parameters of the corresponding sinusoids. A 




(1) 
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transient in the beginning of a segment results in low sinusoidal frequencies, while a transient 
in the end of the segment results m high sinusoidal frequencies. The frequency resolution of 
the sinusoidal model depends on the required resolution in estimation of transient locations. 
If the required time resolution is one sample then the required frequency resolution is defined 
by the reciprocal of the DCT size. 

Due to the duality between the transient location in a time domain segment 
and the frequencies of the corresponding sinusoids, the obvious way to modify the transient 
location is to modify the corresponding frequencies (plus a correction in the phase 
parameters). The transient location in the time domain segment is denoted by no and the 
closest allowed location from a time grid is denoted by n . Then the desired time shift is 
defined as 

An = m-n (2) 

In order to modify the transient location by An the frequencies cOy and phases 
0ij corresponding to the transient should be modified as follows: 

Ann: 

(a,j-m,j . 

N (3) 



^., = ^„+^(i_l + 0-l)f) (4) 
No modification of amplitudes Ay is needed. 

Note that the above procedure is different from independent quantization of 
sinusoidal parameters. All frequencies corresponding to one transient are modified by the 
same amount. This, together with the phase correction of equation (4) above, ensures that the 
shape of the time domain transient is preserved, only the location is modified. 

Because the DCT size is relatively large at one second, more than one 
transient can occur in a time domain segment. In this case, the model has to identify 
sinusoidal parameters corresponding to different transients. This is done by declaring close 
sinusoidal frequencies cOy to represent the same tiransient. Specifically, two sinusoids having 
frequencies differing by not more than e^o are declared to represent the same transient and two 
sinusoids having frequencies differing by more than are declared to represent different 
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transients. Then locations of all transients are modified separately. Below when reference is 
made to a group of frequencies co/j reference is being made to frequencies corresponding to a 
particular transient. 

A transient can occxir at the beginning or at the end of a time domain segment. 
In this case, the modification of sinusoidal frequencies can yield frequencies below 0 or 
above n. This results in the distortion of the shape of the time domain transient. To account 
for this, an overlap is allowed between time domain segments (0.1 seconds). In this case a 
transient can appear in two overlapping segments, i.e. in the region of mutual overlap. 
Because the overlap is sufficiently large, if the transient is located very close to a border of 
one of the overlapping segments, then it is located at a safe distance from a border of the 
other segment. It is straightforward to identify the transient location from sinusoidal 
frequencies, and therefore it is easy knowing the estimated sinusoidal frequencies in the two 
overlapping segments to identify when a transient is represented in two segments. If such a 
situation occurs, the corresponding sinusoids in the segment are cancelled where the 
transient is closer to the corresponding border. 

A typical transient lasts for more than one time sample. A natural question is 
then what is the location of of the transient. After the modification of location the 
corresponding sample of the transient will be placed at location n corresponding to the 
beginning of a segment defmed by the time grid. Therefore, it is important that the estimated 
value no corresponds to the start of the transient. The time domain approach described below 
has proved to yield good results. First, the time samples rimm and rtrmx are identified 
corresponding to the frequency values min(coij) and max((Oij), where coij are frequencies of 
sinusoids corresponding to a particular transient. Next, the highest amplitude of the 
estimated transient signal in the time interval [riram, "max] is found. Then, the start sample of 
the transient no is defined to be the first sample in the interval [«min, «max] having amplitude 
higher than 10% of the highest amplitude. 

Typically, the estimated transient component of an audio signal contains 
samples of small amplitudes before the sample no. Because the time sample no is declared to 
be the first sample of the transient and that no transient can occur at a distance defmed by em 
before the transient, the corresponding samples before no are forced to have zero amplitude. 
As a result, those samples go to the residual signal with their original amplitudes. 

Having estimated the location of transients and modifying their location as 
described above the modified signal can now be modeled to allow the signal to be coded. 
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A damped sinusoidal model is used to model the modified signal, which aims 
at approximating a signal s with a sum of sinusoids with exponentially modulated amplitudes, 
i.e. 

2M 

s(n) = ^ B^e"""" cos(v„n + yf„) 

(5) 

= Y^r^pl,n^O,...,K-\ 

where rm,Pm g C. i<r e iV is the segment length. Equation 5 expresses s(n) as the sum of M 
damped (complex) exponentials. The parameter determines the initial phase and 
amplitude, while pm determines the frequency and damping. In order to determine the 
parameters and pm for the M exponentials the matching pursuit algorithm was used, as 
described in "Matching pursuits with time- frequency dictionaries", IEEE Transactions of 
Signal Processing, Volume 41, pp 3397-3415, December 1993. Matching pursuit 
approximates a signal by a finite expansion into elements chosen from a redundant 
dictionary. Let Z) = (gy)^^, be a complete dictionary of unit- norm elements. The matching 

pursuit algorithm is a greedy iterative algorithm which projects a signal s onto the dictionary 
ry element gy that best matches the signal and subtracts this projection to form a residual signal 
^5 to be approximated in the next iteration. Finding the best matching dictionary element 

consists of computing the inner products {s, gy) and selecting the element that maximises the 
f y inner product. In order to find the parameters and p„ a dictionary is constructed consisting 
ill of damped exponentials, 

g^^ =ce°"e''",n = 0,...,i^-l (6) 
20 Where the constant c is introduced for having unit-norm dictionary elements, 

and compute the inner products of the residual signal at iteration w, Sm and the dictionary 
elements defined in equation 6: 

\K ^Ea,)} = cYa inY"e' , (7) 

By doing this for different values of a, the transfer fimction Sm(z) is evaluated 
25 on circles in the complex z-plane having radius e" 

The method described above has been experimentally tested and the following 
gives results and discussion of computer simulations and informal listening tests performed 
on audio signals. The audio excerpts used were a castanet signal, songs by ABBA, Celine 
Dion, Metallica and a vocal by Suzanne Vega. The signals were sampled at 44.1 kHz. The 



ill 
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DCT size is 44288 samples (approximately 1 second) and the overlap between time domain 
segments is 4410 samples (0.1 seconds). The sinusoidal analysis of the DCT domain signals 
is done using Manning windows of length 512 samples and mutual overlap of 256 samples. 
The transient component of the signal was estimated and subtracted to form the residual 
signal. Next, the transient locations were modified according to a time grid of 220 samples 
(approximately 5 ms). 

It is important to verify that the modification of the transient locations does not 
introduce any audible distortion. To check that, the modified transient signal was added to 
the residual signal. The listening tests conducted verified that there is no perceptual 
difference between the thus obtained signal and the original audio signal. 

In the following, the improvement due to the modification procedure will be 
illustrated. Also discussed is the performance of a damped sinusoidal model with the 
restricted segmentation for an original transient signal (i.e. generally a transient starts at an 
arbitrary location) and the modified transient signal (a transient starts in the beginning of a 
segment). The optimal restricted time segmentation (with the minimum segment size of 220 
samples) for damped sinusoids is found using the technique proposed in "Flexible tree- 
structured signal expansions using time- varying wavelet packets" in IEEE Transactions of 
Signal Processing, Volume 45, pp 333-345, February 1997. The performance is studied in 
terms of signal-to-noise ratio (SNR) versus number of damped sinusoids NDS and is well 
illustrated by Figure 1 where results are presented for a particular transient of the Castanet 
signal; A represents the original transient and B represents the shifted transient. The 
modification procedure results in a considerably smaller number of damped sinusoids needed 
to represent the transient with a certain quality than would previously have been the case. 
Lower plots of Figures 2 and 3, show the reconstruction with 25 damped sinusoids of the 
original and the modified transients, respectively. In these Figures t[ms] denotes time in 
milli-seconds. The original transient is not located m the beginning of the segment and, as a 
result, the modeling error is distributed to samples before the transient. This results in an 
audible pre-echo. On the other hand, the modified transient is located in the beginning of the 
segment and, as a result, the pre-echo problem is eliminated. 

Figure 4 shows a flow diagram of the first embodiment having steps SI to S6, 

where: 

SI represents: Estimate the location of transients in a first time segment of an input signal, by 
a transformation into the fi-equency domain. 



PHNL010529 



12 23.07.2001 

52 represents: Modify the location of the transients in the spatial domain by modifying the 

corresponding frequencies, to locations on a predetermined time scale. 

53 represents: Estimate the location of transients in second and subsequent time segments of 

the transient signal, by a transformation into the frequency domain. 

54 represents: Modify the location of the transients in the spatial domain by modifying the 

corresponding frequencies, to locations on a predetermined time scale. 

55 represents: Decompose an audio signal into transient, tonal and noise components. 

56 represents: Recombine the decomposed signal for transmission or playback. 

It may be possible that a similar improvement to that mentioned above would 
be achieved in the case of a full- precision variable segmentation (and no signal 
modification). However, the restricted segmentation and the modification procedure result in 
a much lower total computational cost. Also, less side information is required to describe the 
restricted segmentation. 

A second embodiment of coding method involves a different method of 
estimating the location of transients in an input signal and a different modification procedure. 
The locations of transients are modified in such a way that a transient can only occur at the 
beginning of a sinusoidal segment, which sinusoidal segments are defined by a specified 
segment size, which may be 5 milliseconds (ms); this is referred to as a restricted 
segmentation, and corresponds to that of the first embodiment. The reference to a beginning 
of a sinusoidal segment can be taken to be a reference to a beginning of a time grid in the &st 
embodiment; the reference to a sinusoid simply refers to the modeling procedure used. 

This second embodiment uses the same idea as the first embodiment in that 
transient locations are modified to improve the modeling of signals, m particular, audio 
signals. However, this second embodiment provides an improved method of modifying the 
location of transients. 

To summarize the first method, the input signal was modified by estimating 
the location of transient components using a model based on the duality between the time and 
jS-equency domain for the signal; subtracting the transient component; modifying the 
locations of transients such that their beginnings can only occur at the beginnings of 
sinusoidal segments and a restricted segmentation; and adding the modified transient to the 
residual signal in order to obtain a modified audio signal. 

In outline, the method of the second embodiment involves detecting the 
beginnings and ends of transient and audio signal using an energy based approach with two 
sliding rectangular windows, as described in "Audio subband coding with improved 
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representation of transient signal segments", from proceedings of EUSIPCO, pages 2345- 
2348, Greece 1998, incorporated herein by reference; followed by moving the identified 
transients to locations specified by a chosen time grid or sinusoidal segmentation grid; and 
time-warping parts of the signal between the identified transients in order to fill the intervals 
5 between the modified transients. 

The transient detection approach as described in "Audio subband coding with 
improved representation of transient signal segments" mentioned above, is based on the 
evaluation of the criterion fimction, C{n) : 



Q 
Q 

m 
m 

\jl where « is a time sample, (n) and E,^ (n) are the energies of the input signal within 
length-7\r rectangular windows on the left- and right-hand side of the time sample n. 

Significant peaks of the criterion function C(n) correspond to the beginnings of transients. 
C3 The end of a transient is defined by searching the first value of C(n) after the beginning of a 

ry 

Pi transient, which is just below a certain threshold. 

Once the beginnings and ends of the transients have been located using the 
20 above method the transients are simply removed firom the signal and relocated to the nearest 

location on the specified sinusoidal segmentation grid, effectively by a cut and paste method. 

This part of the procedure is particularly straightforward and is easily implemented by the 

person skilled in the art. 

As would be appreciated, due to the modification of the transient locations, the 
25 distance between two consecutive transients in an audio signal can become longer (e.g. if one 

is shifted forward and the other is shifted backward), or the distance can become shorter (e.g. 

if a first transient is shifted backwards and a second transient is shifted forwards in time). In 

figure 5 examples of transient modification where the distance is increased is shown, whereas 

in Figure 6, a reduced distance between transients is shown. In order to fill the interval 
30 between the modified transients the signal part in between must be modified in some way to 

allow for the greater or smaller distance between transients. 
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The signal is modified by time-warping, this is done in such a way that 
preserves the correct amplitudes of the edge points of the signal in between the transients, 
thus there are no discontinuities introduced just before or just after a transient, as described 
below. The time-warping results in the signal between transients being stretched (as shown 
in Figure 5) or compressed (as shown in Figure 6). To compute the amplitudes at the new 
integer sampling positions based on the known amplitudes of the original samples, a band 
limited interpolation method based on sine functions is used (the bandlimited interpolation is 
described in Proakis and Manol^s "Digital Signal Processing. Principles, Algorithms and 
Applications", Prentice-Hall Intemational, 1996). Modified Harming window is used. To 
compute the amplitude of each new sample, amplitudes of eight original samples are used, 
four at each side of the new sample. 

The stretching or compressing of a signal results for tonal signals in a 
corresponding change of the fundamental frequency,/o. The goal of the modification 
procedure is to ensure that the induced modifications of fo are not audible. 

In order to achieve the modification, the following algorithm is used for time- 
warping the part of the signal between the two identified and modified transients; 

(a) if the required change in length of a signal part in between two transients results in a 
change of fohy no more thmi 0.2%, the signal is simply subjected to a band limited 
interpolation method based on sine functions. This is the example shown in Figures 5a 
and 6a. If /o changes by more than 0.2% then follow step b) as described below. 

The reason for the limit of 0.2% is that it has been determined from the 
literature on psycho-acoustics that cheaiging fo of a tonal sound by 0.2% can be audible, as 
described in "An introduction to the psychology of hearing", Academic Press, 1997. Our 
own experiments verify this result. 

(b) The signal part is split in between two transients into two non-overlapping intervals; the 
first interval is located directly after the end of the first transient and lasts 10 ms (as 
illustrated by interval 1 in figures 5b and 6b), and the second interval is the remaining 
part, i.e. it lasts until the beginning of the second transient (as shown by interval 2 in 
figures 5b and 6b). The lengths of the two intervals are modified by a different amount. 
If the required change in length of the signal part in between two transients can be done 
by changingy^ in the first interval by no more than 2% and in the second interval by no 
more than 0.2%, then the signal in the two intervals is time-warped correspondingly as 
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shown in the lower parts of figures 5b and 6b. Otherwise go to step c) as described 
below. 

The reasoning behind step b) is that the interval directly after the end of a 
transient is the interval where the masking effect from the tr^sient is strong. Therefore, 
larger changes of the signal in this interval are possible before they become audible. Our 
experiments verify that a change of/oby no more than 2% in the interval 10 ms directly after 
the end of a transient is inaudible. 

(c) time-warp the signal in the two intervals such that the resulting change of /q is no more 
than 2 % in the interval 1 and no more than 0.2 % in the interval 2. If the resulting change 
in length is not sufficient to fill the distance between the shifted transients then apply an 
overlap-add procedure with a modified Banning window using samples from the two 
intervals in order to increase or decrease the length of the signal. To ensure a smooth 
transition between two intervals, the length of the overlap-add region is chosen to be 
larger than required to obtain a correct length of the signal in between two transients 
(figures 5c and 6c). 

In figures 5 and 6 the new locations of transient beginnings are depicted with 
small arrows. In figure 5 the signal part in between two transients becomes longer. In figure 
6 the signal part in between two transients becomes shorter. In the lower part of figure 6c a 
small vertical shift is shown for clarity's sake. 

Various computer simulations of the method of the second embodiment, 
together with informal listening tests with audio signals were carried out. The audio excerpts 
used were castanets, bass, trumpet, Celine Dion, Metallica, harpsichord, Eddie Rabbit, 
Stravinsky and Orff. The signals were sampled at 44.1 kHz. The transient locations were 
modified according to a time grid of 220 samples (approximately 5 ms). It is important to 
verify that the modification of transient locations does not introduce any audible distortion. 
The listening tests conducted verified that there is no perceptual difference between the 
original and modified audio signals. 

Next, it was demonstrated that there is an improvement in the modeling of the 
signal due to the modification procedure. A comparison was made between the performance 
of a damped sinusoidal model with the restricted segmentation for an original transient signal 
(i.e. generally transient starts at an arbitrary location) and for a modified transient signal (a 
transient starts at the beginning of a segment, as defined by the present method). The lower 
parts of figures 7 and 8 show the reconstruction with 25 damped sinusoids of the original and 
the modified transients, respectively. The original trmisient is not located at the begkming of 
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the segment, and as a result, the modeling error is distributed to samples before the transient. 
This results in an audible pre-echo, shown by the amplitude of the signal and the lower part 
of Figure 7 between 5 ms and approximately 7.5 ms, which is not shown in the upper part of 
the Figure 7 that shows the original transient. On the other hand, the modified tramient is 
5 located at the beginning of the segment and, as a result, the pre-echo is eliminated as 

demonstrated in Figure 8 in that the amplitude of the signal for upper and lower parts of the 
figure moves from zero inmiediately after 5 ms, i.e. both at the same time. 

Figure 9 shows a flow diagram of the second embodiment having steps Tl to 

T6, where: 

10 T 1 represents: Estimate the location of transients (beginning and end) in a first time segment 

of an input signal, by an energy based approach. 
iA T2 represents: Modify the location of the transients by cutting and pasting to locations on a 

11 predetermined time scale, and timewarp the signal parts in between. 
S3 T3 represents: Estimate the location of transients (begmning and end) in second and 
^ subsequent time segments of the input signal. 

;p T4 represents: Modify the location of the transients as above, and timewarp the signal parts in 
i! between. 

hi 

T5 represents: Decompose the audio signal into transient, tonal and noise components. 
Jfi T6 represents: Recombine the decomposed signal for transmission or playback. 

The method described in the second embodiment provides a more general 
procedure and provides good results, which are an improvement on those of the first 
embodiment. The time-warping principal is based on the knowledge of sound perception and 
the procedure of the second embodiment is less complex to implement and utilize. 

The advantages of the second embodiment over prior art methods and also the 
first embodiment are that the transient detection model is more general and provides good 
results for various transients, not just short transients. Also, the time-warping of the signal 
parts between transients is based on the knowledge of the properties of sound perception, 
such as pitch perception and temporal masking effects. Furthermore, the method of the 
second embodiment results in a significantly lower computational complexity. 

Both of the methods disclosed herem provide a particularly advantageous 
method for coding audio and video signals. In particular, restricting the transient locations 
simplifies the analysis procedure in an audio coder (involving transient, sinusoidal and noise 
models) significantly. Also, the side information associated with the corresponding 
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segmentation is reduced because of the restricted segmentation often used in the two 
embodiments described. 

Furthermore, the introduced difference in transient locations is not of 
perceptual importance. 

The method could be implemented in devices for storing, transmitting, 
receiving, or reproducing audio and/or video, e.g. solid state audio devices. Figure 10 shows 
an audio coder 10 and an audio decoder 12 which receive an audio si^al (A) for coding and 
a coded signal (C) for decoding respectively, with the decoder 12 outputting the audio signal 
A. In particular, the audio coder may be included in a transmitting or recording device, 
further comprising a source or receiver for obtaining the audio signal and an output unit for 
transmitting/outputting the coded signal to a transmission medium or a storage medium (e.g. 
a sold state memory). For stereo audio signals, the time and intensity with which a signal 
reaches both ears play a major role on localization of sounds, i.e. the perception of direction 
and distance to the sound source. More precisely, it is the difference in time (interaural time 
difference) and difference in intensity (interaural intensity difference) with which the signal 
reaches both ears, which form the so called stereo image. Here, we deal with time 
modifications of audio signals for the purpose of efficient modeling. Therefore, below we 
will concentrate our attention on the resulting interaural (interchannel) time differences. 

The audibility of interchannel time difference and relative importance of 
transients and ongoing parts in formation of stereo image depend upon a variety of factors, 
including duration of sounds, frequency content, repetition rate (for transients). The 
important result, however, is that interchannel time differences as small as of order of 10 us 
can be detected by the auditory system (using cues either from transients or ongoing parts). 

When modifying transient locations, also the ongoing parts are modified due 
to the time shift and time warping, i.e. both important cues are present. Therefore, care has to 
be taken for not destroying the original stereo image. 

An efficient modeling with damped sinusoids can be obtained if transient 
locations in both stereo channels are modified such that the transients start at the beginnings 
of the sinusoidal segments. The independent modifications in the two channels would, 
however, generally result in a destroyed stereo image. A possible solution to this problem 
could be to modify the transient locations according to the sinusoidal segmentation before 
modeling with damped sinusoids, but to send side information describing the original time 
differences between corresponding transients in the two channels to the decoder. The , at the 
decoder the synthesized signal in one of the channels can be unwarped according to the 
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original time difference. As a result, the synthesized transients occur generally at locations 
different from their original locations but the interchannel time difference between the two 
transients is preserved. This solution is especially suitable for highly-correlated stereo 
channels, having similar detected transients with low interchannel time differences. 
5 It should be noted that the above-mentioned embodiments illustrate rather than 

limit the invention, and that those skilled in the art will be able to design many altemative 
embodiments without departing from the scope of the appended claims. In the claims, any 
reference signs placed between parentheses shall not be construed as limiting the claim. The 
word 'comprising' does not exclude the presence of other elements or steps than those listed 

10 in a claim. The invention can be implemented by means of hardware comprising several 
distinct elements, and by means of a suitably programmed computer, hi a device claim 
enumerating several means, several of these means can be embodied by one and the same 

K item of hardware. The mere fact that certain measures are recited in mutually different 

Q dependent claims does not indicate that a combination of these measures cannot be used to 

Q 

advantage. 

^3 In summary, an improved representation of trmisients in audio signals 

ry comprises modifying transient locations in such a way that a transient can occur only at a 

beginning of a sinusoidal segment. The modification procedure comprises the steps: 
!f : - detecting a beginning and an end of a transient using an energy-based approach with two 
two sliding rectangular windows; 

- moving samples between the beginning and the end of the transient to the locations 
specified by the segmentation used; and 

- time- warping the signal parts in between the transients in order to fill the intervals between 
the modified transients. 



